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' Abstract: It is shown that the flavour projection of staggered fermions can be written 

^ ' as a projection between the fields on four separate, but parallel, lattices, where the fields ^ 

psj ■ on each are modified forms of the standard staggered fermion field. Because the staggered ^ — ^ 
\^ • Dirac operator acts equally on each lattice, it respects this flavour projection. We show that 
the system can be gauged in the usual fashion and that this does not interfere with flavour 
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Q , projection. We also consider the path integral, showing that, prior to flavour projection, it 

evaluates to the same form on each lattice and that this form is equal to that used in the 
quarter-root trick. The flavour projection leaves a path integral for a single flavour of field 



$H ' on each lattice. 
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1. Introduction 

Constructing lattice field theories that incorporate a single non-degenerate fermion field, 
with chiral symmetry, has proven to be very challenging. Wilson originally identified 
the problem by showing that naively discretizing the continuum Dirac equation led to a 
doubling in the number of fermion fields in each dimension of momentum space [||. As 
a solution, he proposed adding a term proportional to the second derivative to break the 
degeneracy between the physical field and the doublers. However, this term broke the chiral 
symmetry of the theory. Subsequently, Kogut, Susskind and Banks proposed the staggered 
fermion formulation QlQQ in which the number of degenerate fields is reduced to 2"/^ 
in n dimensions and which retains a limited form of chiral symmetry. Other formulations 
have since succeeded in isolating a single chiral field by introducing an infinite number of 
regulating fields In the finite case, these formulations contain an exponentially 

suppressed term that breaks the chiral symmetry of the theory. 

For reasons of computational efficiency, the staggered fermion formulation has re- 
mained popular in lattice calculations and it is common to see the quarter-root trick used 
to model a system with one flavour of fermion Q j^. It has been shown that, in the free field 
case, it is possible to construct a lattice Dirac operator whose determinant is the quarter- 



root of the determinant of the conventional staggered Dirac operator |1C]. Although, the 
locality of this operator has been questioned [0], numerical studies suggest that a suitably 
local operator does exist [^]. 

In this paper, we draw on previous work [^] to show that the flavour projection of 
staggered fermions can be written as a projection between four modified staggered fermion 
fields on four separate, but parallel, lattices. The staggered Dirac operator respects the 



flavour projection because it acts upon the fields on each lattice equally. We show that 
we can gauge the fields in each lattice in the usual way and that the flavour projection 
is equally valid in the gauged case. We discuss the role of chiral symmetry and we see 
that, in this framework, the 75 (8) 75 chiral symmetry of the staggered fermion formulation 
can be used to isolate the chiral components of the fields. Finally, we consider the path 
integrals on each lattice, before and after flavoTir projection, showing that projection leaves 
the path integral for a single flavour of field on each lattice. In considering this, we sec that 
evaluating the path integral, prior to fiavour projection, leads to an expression identical to 
that used in the quarter-root trick. 



2. Background 

Our starting point is the staggered fermion formulation in four Euclidean dimensions. The 
lattice Dirac equation is diagonalised by rewriting the spinor field as ipn, where 

and 



massless action is 



We retain only the first component of -^n, defining the field Xn = V'n- '^^^ resulting 



where 



2 

n,n 



r/i(n) = 1 

ry^(n) = : for /i 7^ 1 . 

The lattice sites are grouped into hypercubes of side twice the original lattice spacing, 
each containing sixteen sites. We relabel the coordinate of the lattice site n with the 
coordinate of the origin of the hypercube, m, and the location of the site within the 
hypercube, s. The relationship between the four component vectors n, m and sisn = m+s. 
Together the hypercubes yield a fermion field that is free from degeneracy in momentum 
space 

^m=Afo^T{s)xtn, (2.1) 



where A/q is a normalisation constant. For TVq = the action is 



(7M^^)^M + ^(75^7;^75)nM 



i>m , (2.2) 



m,fj, 

where A^B is the product of the spin and fiavour spaces and b is twice the original lattice 
spacing. 

^pm is now a four by four matrix containing four degenerate flavours. The action mixes 
both the spin and flavour space and the only operator that anticommutes with both terms 
in the action is 75 (8> 75, generating a limited form of chiral symmetry. 



3. Flavour projection 

In the continuum limit, the second term in equation (|2.2| ) goes to zero, we have that 
V'm ip{m) and the columns of il){m) come to represent separate flavours of the field. 
Each column can be isolated with the operator P^''), where 

and the projection matrix P*-^^ is of the form 

diag(l, 0, 0, 0) : for 6 = 1 , diag(0, 1, 0,0) : for 6 = 2 , 
diag(0, 0, 1,0) : for 6 = 3 , diag(0, 0, 0,1) : for 6 = 4 . 

P^*^ is defined through 

p(^) = i (1 + iab7i72) (1 + /3671727374) (3.1) 

with 

Ofe = (-1,+1,-1,+1)^ 
/3b = (-1,-1,+1,+1)^ . 

We can expand the projection matrix to give 

p(^) = ^ (1 + iafe7i72 + /5fe7i 727374 - iabPhlsn) (3.2) 

and this is equal to 

P(^) = i (diag(l, 1, 1, 1) + abdiag(-l, 1, -1, 1) 

+;a6diag(-l, -1, 1, 1) + a6/3ftdiag(l, -1, -1, 1)) . 

Writing 

Pi = diag(l, 1, 1, 1) , P2 = diag(-l, 1, -1, 1) , 

Pa = diag(-l, -1, 1, 1) , P4 = diag(l, -1, -1, 1) , 

we can see that [P',P*^''''] = 0, where D is the staggered Dirac operator, because 

[Dij{m)]P, = D[ij{m)Pi] . 

At finite lattice spacing, 6^0 and the second term in the staggered Dirac operator 
breaks this associativity. Thus [-DjP] 7^ 0. However, if we could transform the Pi in a 
manner consistent with the effect of D on ipm , we would be able to maintain [D, P^^^] = 
at finite lattice spacing. 

To implement this strategy, we write ■4^mP^^^ as 

V'„^P(^) = l^Pn, (Pi + abP2 + PbPs + abPbP4) 

= i (^m + "b^m + Pb'^fa + abPb'^t) ■ 

We introduce three additional lattices, in parallel to the first, and on the i-th lattice, Li, 
we place the field ^1^. Flavour projection now becomes a projection between the lattices 



and to describe it we introduce the following operator to map a field from lattice Lj to 
lattice Lj 

Writing the fermion field across all four lattices as = X^j using equation (^.Ip, 

the flavour projection operator can now be written as 

PW^^ = I (l + /3b (T31 - Tia + T42 - T2a) ) 

x (l + ab (T21 - T12 + T43 - T^a) ) ■ 

Application of p(^)^' will leave column i on Lj. However, because D mixes the spin 
and flavour spaces equally on all four lattice, the same degrees of freedom of ipm will be 
projected out, irrespective of the order in which we apply D and P^^^^ to i/jm- Thus we have 
[PW,Z)] = 0. 



4. Chiral symmetry 

With the components from only one flavour of field on each lattice, we can see that the 
residual 75 (8)75 chiral generator is sufficient to isolate the chiral components. For example, 
the application PfiP^^^^ leaves the positive chirality components of the first column of 
1pm on Li and of the second column of ipm on L2. It also leaves the negative chirality 
components of the third column of ipm on L3 and of the fourth column of tprn on L4. 
Because the chiral projection isolates the same components of on each lattice, it does 
not interfere with the flavour projection and so we have that [Pj^/^, P^'')] = 0. 

The effect of the projection on the chiral anomaly can be understood by considering 
P'^^^Pr^. The application of Pr to projects out two positive and two negative chirality 
spinors on each lattice which ensure that the chiral current is conserved. Writing P^^^ as 

P(^) = P^'W 



p^') = i(i + /3,(r3i-ri3 + r. 
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pi*^ = i(i + ab(r2i -ri2 + T43 

we can see that Pa ^ removes one positive and one negative chirality spinor from each lattice, 
leaving one positive and one negative chirality spinor behind. At this stage the chiral 
current is still conserved on each lattice. The subsequent application of P^^"*, separates 
the remaining spinors, placing them on different lattices. This results in a chiral current 
flowing between Li and L3 and between L2 and L4 which manifests itself as an anomaly 
on each lattice. 



5. Gauging the links 

The gauging of staggered fermion fields takes place in the Xn basis. However, the modified 
form of tpm that we place on each lattice leads us to introduce modified forms of xfn- 



Using the identity rr(T(n)tr(m)) = 45"''^, the inverse of equation (|2.1| ) can be shown 
to be 

We define the modified field for lattice i as 

< = ^Tr (r(.)t^g . (5.1) 

By writing the diagonal matrices in terms of the defining 7-matrix combinations from 
equation (|3.2| ), we can relate to Xm follows. 



^Tr {T{s)^m) = Xfn 

'^^^ = w^Tr (T(s)tv;„7i727374) = P(s,i234)X^^^^*' 

'^m = "4^^^ (r(s)"^^m7374) = ip(s,34)Xm'"' • 

Here P(^m,n) is defined to be (—1)*^, where i' is the number of pairs (m,n) with m G M, 
n & N and m > n [^]. We have also introduced C as a complementarity operator. C12S 
gives a vector that complements s in the {1, 2} subspace, but is the same as s in the {3, 4} 
subspace. Similarly, C34S gives a vector that complements s in the {3,4} subspace, but 
is the same in the {1,2} subspace. C1234S gives a vector that complements s in all four 
dimensions. 

On each lattice is gauged in the same manner as Xm- Because of the linearity of 
the transform in equation (5.1), we can apply the projection operator P^^^ to liJ = ^ 



.IS 

ims m 



to project a single gauged spinor onto each lattice. 



6. Path integral formulation 

The path integral contains contributions from all four lattices. On each lattice, the path 
integral takes the form 

[d-^ ][d-^ ]e V y 1 (6.1) 

where Lq is a constant yet to be determined. The fields on each lattice share the same 
degrees of freedom, but, because each lattice lies in a separate space, we can separately 
evaluate the contribution from each. 

The number of elements of ^'^ and that are the negative of their counterparts 
in '^^ is even and this gives the following equivalence 

[di!\d¥] = [di!^\[di!^\ = [d^][di;] , 




giving us the same fermionic measure of integration on each lattice. 



If we consider the action on each lattice written as a matrix equation, it takes the form 

~ i 

5* = "^^KapiA)^^^, where Kafs is the gauged staggered Dirac operator. By rewriting the 
action as 

we introduce iir^^(^) in which an even number of columns are equal to the negative of 
their counterpart in Kap{A). Evaluating 



gives 



{det[K\A)]f' H 

o 



which is equal to {det[K {A)]y\ Because Klp{A), Klp[A) and K^p^A) differ from Klp{A) 
by an even number of columns multiplied by —1, we have that 



(6.2) [SJ 



= {det[K\A)]Y' 
= {det[K{A)])^'> . ^ 

o 

Hence, the products of all four contributions to the path integral is 

{det[K{A)]f'^- . ^ 

o 

For the path integral to maintain the correct continuum limit, we require that Lq = ^. This | 
gives an expression on each lattice that is of the same form as that used in the quarter-root qq 
trick (equation (|6.2D). 

If we evaluate the path integral in the uj basis, we find from equation (|5.2| ) that, for 
each site m on each lattice i, the sixteen components of contain an even number of 
terms that are the negative of their counterpart amongst xfn- Hence, we have that 

[dLU^][du;'] = [dCo'][du;'] = [dx][dx] ■ 

As in the V basis, the —1 factors within the action can be incorporated into the Dirac 
operator to show that the path integral evaluates to {detlM {A)])^° on each lattice, where 
M{A) is the Dirac operator in this basis. 

If, instead of immediately evaluating the path integral on each lattice, we first apply 
P*-*^ to equation ( |6.1D , the flavour projection will affect both the action and fermionic 
measure of integration. As discussed, on each lattice, it will leave the action acting on 
a single flavour of field. It will also project between the measures of integration on each 
lattice, to leave the components belonging to the same single flavour as in the action. The 
resulting path integral on each lattice, describes a form of single gauged fermion. 
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7. Conclusion 



In this paper, we have shown that it is possible to describe flavour projection in staggered 
fermions as a projection between four lattices where modified forms of the fermion field 
occupy each lattice. The staggered Dirac operator respects this flavour projection because 
it mixes the spin and flavour spaces of the fields on each lattice equally. We saw that, after 
flavour projection, the 75 (8> 75 chirality of the staggered fermion formulation is sufficient 
to isolate the chiral components of the remaining fields. We also showed that this system 
could be gauged in the usual way and that it was possible to project out a single gauged 
field. We showed that the projection could be applied to the full path integral and that, 
prior to flavour projection, evaluating the path integral lead to a form consistent with the 
quarter-root trick commonly seen in lattice QCD. 

It is our hope that this work will contribute to the discussion on fermion doubling and 
that, by flnding a natural home for the quarter-root trick, this formulation can contribute 
to the current debate on its validity. 

References 



[1] 
[2] 
[3] 
[4] 
[5] 
[6] 

[7 

[8: 

[9 
[10: 

[11 

[12 
[13 
[14 



K G Wilson, Phys. Rev. D 10 (1974) 2445 



J Kogut and L Susskind, Phys. Rev. D 11 (1975) 395 



T Banks, L Susskind and J Kogut, ^hys. Rev. D 13 (1976) 1043 
L Susskind, \Phys. Rev. D 16 (1977) 3031] 



D B Kaplan, Phys. Lett. B 288 (1992) 342, tiep-lat/9206013 



V Furman and Y Shamir, Nucl. Phys. B 439 (1995) 54, tiep-lat/9405004 



H Neuberger, Phys. Lett. B 417 (1998) 141, hep-lat/9707022 



C T H Davies, \Phys. Rev. Lett. 92 (2004) 02200l| , [hep-lat/0304004 
C Aubin, ^hys. Rev. D 70 (2004) 11450l| , |hep-lat/0407028 



D H Adams, \Phys. Rev. D 72 (2005) 114512| , |hep-lat/0411030 



B Bvmk, M Delia Morte, K Jansen and F Kncchtli, Y^ucl. Phys. B 697 (2004) 343| , 



faep-lat/040302; 



F Maresca and M Peardon, (2004), |hep-lat/0411029| 
S Watterson, PhD Thesis (2007) 



P Becher and H Joos, Z. Physik C 15 (1982) 34S 



